Enhancement of phonon effects in photoexcited states of one-dimensional Mott 

insulators 
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We examine how the electron correlation affects the electron-phonon (EP) interaction in the lin- 
ear optical absorption spectrum of the one-dimensional (ID) extended Hubbard-Holstein model. A 
density matrix renormalization group (DMRG) calculation shows that the effect of the EP inter- 
action on an exciton is enhanced by increasing the on-site Coulomb repulsion. This enhancement 
is in contrast to the effect of the EP interaction on the ground state where the Peierls instability 
is suppressed by the on-site Coulomb repulsion. The DMRG data with the EP interaction fit with 
absorption experiments in ID cuprates better than those for the extended Hubbard model. 
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The interplay of strong electron-electron correlation 
and electron-phonon (EP) interaction has attracted much 
attention [![. One of the issues is to understand the ef- 
fect of the EP interaction on carriers doped into a Mott 
insulator. In the ground state of the Mott insulator, 
the on-site Coulomb repulsion makes the charge distri- 
bution uniform, and then the EP interaction is irrele- 
vant. However, the mobile carriers created in the Mott 
insulator by chemical and/or photo dopings induce the 
charge fluctuation, and thus the EP interaction is ex- 
pected toplay an important role in various physical quan- 
tities H i, II !, S, 0, B i . 

In this paper, we address this issue in the photoexcited 
states of the one-dimensional (ID) Mott insulators. The 
effect of the EP interaction on the photocarriers can be 
observed by optical measurements. In the Raman spec- 
troscopy for Cai.gSro.2Cu03, Raman-forbidden phonon 
peaks are resonantly enhanced for charge-transfer exci- 
tation This observation may indicate the change 
of the selection rule by coupling between phonons and 
photocarriers created by the charge-transfer excitation. 
One of the phonon modes has the local symmetry equal 
to a breathing mode whose effect on the single-particle 
excitation spectra has been intensively studied in high-T c 
cuprates [rH H3|. 

In the photoexcited state of the Mott insulator, there 
are two types of photocarriers called 'holon' and 'dou- 
blon' representing empty and doubly occupied sites, re- 
spectively 13]. Usually, a holon and a doublon form an 
exciton due to the intersite Coulomb repulsion 0, [ll| . 
However, this exciton is different from an electron-hole 
pair in semiconductors, since the strong electron corre- 
lation prohibits the exchange of holon and doublon. It 
has been pointed out that this prohibition enhances a 
magnitude of the third-order optical nonlinear suscepti- 
bility [HI [ItJ ■ Thus, it is important to understand how 
the exciton in the Mott insulators is modified by the EP 
interaction. 

We examine the effect of the EP interaction on the 



exciton in the ID extended Hubbard-Holstein model. 
We perform the density matrix renormalization group 
(DMRG) calculation of the optical spectra. The DMRG 
calculation shows that the effect of the EP interaction 
on the exciton is enhanced by increasing the on-site 
Coulomb repulsion. This enhancement is in contrast 
to the effect of the EP interaction on the ground state 
where the Peierls instability is suppressed by the on-site 
Coulomb repulsion. The DMRG data with the EP in- 
teraction fit with absorption experiments in ID cuprates 
better than those for the extended Hubbard model [li ]. 

The coupling between an electron and a breathing 
phonon in the cuprates can be mapped onto the Holstein- 
and Peierls- type interactions [l9j]. Since the Holstein in- 
teraction is stronger than the Peierls one, we start with 
the ID Hubbard-Holstein model. In addition, we intro- 
duce the nearest neighbor Coulomb repulsion which leads 
to the formation of an exciton. The Hailtonian is defined 
by 
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where cj (cj lCT ) is a creation (annihilation) operator of 

an electron at site i with spin a, and b\ (6,) is a creation 
(annihilation) operator of a phonon at site i. This model 
includes electron hopping, t, on-site and nearest-neighbor 
Coulomb repulsions, U and V , respectively, phonon fre- 
quency, luq, and EP coupling, g. Here, we neglect the 
dispersion of the phonon in order to make our discussion 
simple. 

We calculate the current-current correlation function 
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defined by 



1 



E - H + i-y 



J 



, (2) 



where |0) is the ground state with energy Eq, L is 
the system size, the current operator j is defined by 

j = —itJ2i a( c l o-Ci+i-c ~ H.c), and 7 is a small positive 
number. The optical absorption spectrum is proportional 
to xM/w 2 . 

In order to calculate x( w )i the infinite-system DMRG 
method [2pj is applied to an open boundary system. We 
have also carried out the finite-system DMRG calculation 
for some data, and this calculation quantitatively agrees 
with the infinite-system DMRG data. In eq. fl]), we have 
introduced the y-term, V(rii — l)(rii + i — 1), instead of 
the standard notation, VriiUi+i, and the present notation 
does not produce the boundary effect which is localiza- 
tion of holon and doublon at the edge. We introduce 
a mixed-state density matrix which is composed of the 
ground state |0), a photoexcited state j* |0), and two cor- 
rection vectors (e+E — H |0) with e = ui, uj+2^. 
In terms of the reduced Hilbert space made by the den- 
sity matrix, the spectrum for lo < e < a; + 27 is calculated 
after one DMRG run. The number of the eigenstates of 
this density matrix are taken up to 300 for renormaliza- 
tion. In order to correctly deal with the phononic degrees 
of freedom, we use the pseudo-boson method [2l|. In this 
method, a boson operator with a restricted Hilbert space 
is exactly transformed into a set of hard-core bosons. For 
instance, the boson operator with 4 states is transformed 
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the hard-core boson operators a\ and 0,2- Then, step-by- 
step renormalization of each boson is carried out. For 
g < \/0.5i, we have taken 4 hard-core bosons per site, 
which correspond to 16 phonon states per site. 

In numerical calculations, we take L = 20, 20 electrons 
(half-filling), U = lOt, iv = 0.5i, and 7 = 0.2t < u . 
The magnitude of ujq may be three-times larger than the 
frequency of the breathing mode in cuprates [l9[ in or- 
der to make numerical calculation feasible. However, in 
the Hubbard-Holstein model, a dimensionless parame- 
ter, A = g 2 /AtujQ, is more important for estimation of 
the strength of the EP interaction, and thus A is used 
for comparison of numerical data with observed optical 
spectra. 

In Fig. [TJ we show x( w ) f° r various g and V values. 
For 5 = 0, the V dependence of the spectra has been 
studied intensively, and the present data are consistent 
with the previous data [15[ : for V < 2t only a continuum 
band exists, which is composed of nearly free holon and 
doublon, while for V > 2t the holon and doublon form an 
excitonic bound state. A pronounced peak at lu = 6. It 
in Fig. 1(d) and a tiny hump at u) = 6M in Fig. 1(c) 
correspond to the exciton (the other spiky structures in 
Figs. l(a)-(d) come from the finite-size effect). For g 7^ 0, 
the DMRG data show that the effect of g on the spectra 
becomes remarkable with increasing V, and in particular 




FIG. 1: Current-current correlation function xi 1 - ) f° r the 
Hubbard-Holstein model with U = lOi. (a) V/t = 0, (b) 
V/t = 1, (c) V/t = 2, and (d) V/t = 3. 



the bound exciton for V > 2t is strongly renormalized by 
the EP interaction. 

For V < 2t and j ^ 0, the band edge shifts toward 
lower-energy region. This is because the holon and dou- 
blon are dressed with phonons, leading to an energy gain 
of the photoexcited state. The energy gain estimated 
by the DMRG data is smaller than 2g 2 /u>o given by the 
second-order perturbation expansion with respect to g 
for t = 0. For instance, the energy gain for V = t and 
g = \/0.5t in Fig. 1(b) is about 0.5i, and this is smaller 
than 2g 2 /u>o — 2t. In this case, the motion of holon and 
doublon suppresses the effect of g on these carriers. 

For V > 2t and j / 0, the exciton peak for g = 
splits into multiple peaks. For g — \/0.5t in Fig. HKd), 
these peaks are positioned at 5.38i and 5.85t. The en- 
ergy difference between them is 0.47i ~ ujq, and thus 
it is expected that the exciton has a polaronic feature. 
Because of the splitting, the line shape of the exciton is 
asymmetric (see Fig. |TJd)) . The shape is different from 
that obtained by a classical treatment of phonons [22| . 
The energy gain of the exciton by g is still smaller than 
2g 2 /ojQ, even though the exciton would be localized for 
V = 3t. 

Next, we examine whether the effect of g on the ex- 
citon is enhanced or reduced by U. For this purpose, 
we compare the spectra for U = lOi and in the large-?/ 
limit [23j. It is usually recognized that the parameter 
U = lOt is characteristic of strong electron-electron cor- 
relation. Actually, for g = 0, the absorption spectra for 
these parameters are almost the same. In this case, the 
exciton is bounded for V ^ 2t and U = lOi, and this crit- 
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ical V value is close to the critical one V — 2t in the 
large-?/ limit [HI]. However, the line shape of the stress 
tensor correlation function changes dramatically even for 
lOt < U < 40t [HI. As will be shown later, we can 
see quantitative difference between the absorption spec- 
tra for U = Wt and in the large-?/ limit in the presence 
of g. This means that the effect of g on the spectra is 
sensitive to the value of U. 

In order to study the optical spectra for U — > oo, we 
introduce a holon-doublon model coupled with phonons. 
This model describes the charge sector of the photoex- 
cited states in the large- U limit of the Hubbard-Holstein 
model ([I]). We drop the spin degrees of freedom be- 
cause of the spin-charge separation in the ground state, 
and thus the numerical calculation becomes easier. The 
Hamiltonian is defined by 
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where h\ (hi) and d\ (di) are creation (annihilation) op- 
erators of holon and doublon at site i, respectively. The 
holon and doublon obey the local constraint described by 
h\hi + d\di = 1. The current-current correlation function 
can be described by 



ttL 



-Im 1 



1 



(4) 



where |1) is a state with one holon-doublon pair. This 
state is created by applying the current operator j — 



i+U 



to the classical Neel ordered 



state. Here, we introduce the prefactor 2 in order to 
take account of the effect of quantum spin fluctuation. 
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the lowest-energy photoexcited state is an exciton. In 
order to examine the difference of the phonon effect on 
the exciton between these models, g is changed from to 
y/0.5t. For instance, the difference appears in Fig. [DJc): 
the spectral weight of the lowest-energy peak is largest 
for U = lOt, while the weight of the lowest-energy peak is 
smaller than that of the second-lowest-energy peak in the 
large-?/ limit. Therefore, the increase of U significantly 
changes the lineshape of the exciton in the presence of 
g. This means that the effecive EP interaction increases 
with U. 

Let us consider why U enhances the effect of the EP 
interaction on the exciton. First of all, it would be useful 
for later discussion to introduce the exactly solvable case 
of the holon-doublon model. In the large- limit, by the 
canonical transformation anologous to the method taken 
in ref [24j , we have the following result 
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where Ai = 87r(2g 2 /w§) ; exp (~2g 2 /u>%)/ll and ^ = 
U -V - U 2 /V - 2g 2 /u + luj . We find that the spec- 
tral weight for the lowest-energy photoexcited state |no), 
Aq = 87rexp (— 2g 2 /uj 2 ), decreases exponentially as g in- 
creases. On the other hand, the decrease of the spectral 
weights for higher-energy photoexcited states \ni) is sup- 
pressed by a factor (2g 2 /w 2 ) 1 . When g > ujo/V2, Aq 
becomes smaller than A\ because of strong renormaliza- 
tion of |no) due to large EP coupling. In this case, we 
suppose large expectation value of the phonon number 
for |no), Nq — (no\ J2i b\bi \no). In the following consid- 
eration, we assume this parameter region where |no) is 
strongly renormalized and A < A±. 

In the large-{7 limit, charge fluctuation does not occur 
in the ground state, and |no) is composed of the bases 
with one holon-doublon pair. Because of the inhomo- 
geneous charge distribution, phonons are spontaneously 
excited around the pair, leading to large A^o- On the 
other hand, when U is finite, the bases without the pair 
exist in |no). Since these bases represent homonegeous 
charge distribution, N decreases as U decreases. 

Next, we examine Aq as a function of No in order to 
understand the relation between Aq and U. It is clear 
that the dominant contribution to the ground state comes 
from |0) ~ \AF) |0), where \AF) and |0) represent the an- 
tiferromagnetic state and phonon vacuum, respectively, 
as long as the ground state is the Mott insulator. In this 



case, Aq 



\no\J 



is given by 



FIG. 2: Comparison of x( w ) between the Hubbard-Holstein 
model with U = Wt (fine solid line) and the holon-doublon 
model (bold solid line) . We take V = 2t. The spectra for the 
holon-doublon model are shifted for the comparison. 

Figure [2] shows comparison of x( w ) between the 
Hubbard-Holstein model with U — lOt and the holon- 
doublon model. Here, V is taken to be 2t, and then 
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where the holon-doublon pair |0, fl) is located on the i- 
th and (i + l)-th sites. This equation shows the relation 
between Aq and Nq. When Nq decreases, the weights 
of the bases without phonons become dominant in \uq), 
and this effect increases Aq. According to the previous 
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paragraph where we found that Nq decreases with de- 
creasing U, A increases with decreasing U. We have 
also confirmed U dependence of No and A$ by exact di- 
agonalization on small clusters, and the diagonalization 
result agrees well with the above consideration. 

Therefore, it is concluded that the effect of the EP in- 
teraction on an exciton is very strong in correlated elec- 
tron systems. This is essentially different from the ef- 
fect of the EP interaction on the ground state where the 
Peierls instability is suppressed with the increase of the 
on-site Coulomb repulsion. 

Here, we compare the present DMRG results with 
the optical absorption spectrum in the ID cuprate 
Sr2Cu03 [HI]. The absorption spectrum has asymmet- 
ric structure: the high-energy side of the spectrum is 
broad, while the low-energy side is fitted with a single 
Lorentzian function with a broadening factor of the or- 
der of 0.1 eV. In addition, we find a flat part at the top of 
the absorption peak (1.75eV $ ui 1.82eV). The energy 
range 0.07 eV of this flat part is comparable to a phonon 
frequency of the breathing mode seen in Cai.sSro^CuOs. 
We consider that the flat part is a signature of multiple 
phonon peaks, although the clear peak structures would 
be washed away by the phonon dispersion. In our pre- 
vious DMRG calculation for the single-particle excita- 
tion spectra, the dimensionless coupling parameter is es- 
timated to be A - 0.25 (uj = 0M and g = y^t) [|. 



Figure [TJc) or [IJd) with this parameter A is consistent 
with the absorption spectrum. It is noted that V may 
be larger than that estimated by the model without g, 
since the EP interaction makes the spectrum broad. This 
fact might affect the nonlinear optical properties in the 
cuprate [Til ]. 

In summary, we have examined the effect of the EP 
interaction on the exciton in the ID Mott insulator. We 
have performed the DMRG calculation of the optical 
spectra in the ID Hubbard-Holstein model at half filling. 
We have found that the effect of the EP interaction on the 
exciton is enhanced by increasing the on-site Coulomb re- 
pulsion. This enhancement is in contrast to the effect of 
the EP interaction on the ground state where the Peierls 
instability is suppressed by the on-site Coulomb repul- 
sion. The DMRG result with g is consistent with linear 
absorption spectrum of the ID cuprate Si'2Cu03. 
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